Dissipative dynamics of an extended magnetic nanostructure : Spin necklace in a 

metallic environment 
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We study theoretically the dynamics of an "xxz" spin necklace coupled to a conduction electron 
sea, a model system for a nanostructure in a dissipative environment. We extract the long-time 
behavior via a mapping to a multichannel Coulomb gas problem followed by a scaling analysis. The 
strong quantum fluctuations of the necklace cause a nontrivial dependence of couplings on system 
size which we extract via an analysis involving the "boundary condition changing operator", and 
confirm via a detailed numerical evaluation of one case. 



The dissipative dynamics of a magnetic nanostructure 
is important from both technological and fundamental 
points of view. For example, magnetic recording involves 
the polarization of small domains, whose stability over 
long times is crucial. Tunnelling of isolated magnetic 
particles has been extensively studied [1] [2] . We shall be 
interested here in the new physics brought by coupling to 
a conduction electron bath. There are two effects. Ex- 
change of spin between the conduction electron system 
and the nanostructure changes the magnetic state of the 
nanostructure, leading effectively to a tunneling process. 
On the other hand, dissipation from the particle-hole ex- 
citations of the metal suppresses tunnelling [3]. A finite 
cluster of spins coupled to a conduction electron bath is 
thus a model system for examining general issues of quan- 
tum coherence and dissipation in a system with many 
degrees of freedom. Also, recent experimental advances 
in atom manipulation [4] [5] mean that clusters of atoms 
(eg. quantum corrals) can now be placed in controllable 
arrays on suitably chosen solid surfaces, so it seems likely 
that in the near future controlled magnetic nano-arrays 
may be constructed and studied. 

In the case of a single spin in a metal, the inter- 
play between tunnelling and dissipation gives rise to the 
logarithmic renormalizations characteristic of the one- 
impurity Kondo effect [6]. While the two [7] [8] and 
three [9] impurity Kondo problems have been studied, 
there is less known about the behavior of systems con- 
taining larger numbers of spins. The subject has been 
recently studied in the context of the "Griffiths phase" 
scenario for the non-Fermi-liquid physics apparently ob- 
served in certain 'heavy fermion' metals [10]. To study 
this physics, Castro-Neto and Jones [8] considered a clus- 
ter of N magnetic impurities, coupled strongly to each 
other and weakly to a conduction electron bath. They 
presented a series of qualitative arguments and approxi- 
mate mappings from which they concluded that the sys- 
tem could undergo a more or less conventional Kondo 
screening (which they termed a "cluster Kondo effect"), 
at an N— dependent temperature which they estimated. 
The qualitative nature of the treatment suggests how- 
ever that further analysis would be useful. One of us, 



with Morr and Schmalian [11], treated a similar problem 
of N coupled moments from a functional integral point of 
view. In the Ising case, orthogonality effects were found 
to suppress tunnelling in all but the smallest nanostruc- 
tures, while no definitive results were obtained in the 
continuous symmetry case. 

In this paper we attack the problem from a different 
point of view. We identify a specific nanostructure, the 
finite-size spin necklace depicted in Fig.l. There are three 
motivations for this choice of cluster. First, it is (as can 
well be seen) amenable to a detailed analysis. Second, 
it exhibits large (but tractable) quantum fluctuations 
whereas in other works [8] [11] the quantum fluctuations 
within the cluster are not taken into account. Third, we 
can treat the continuous symmetry case as well. 




Fig.l Sketch of 'spin necklace' in conduction electron bath for 
N = 5. Heavy dots represent local moments, which are coupled 
via an exchange coupling K. Shaded area represents conduction 
electron (c) sea, Kondo-coupled to local moments via exchange J. 



The Hamiltonian describing the physics of interest is 
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with Rj the position of spin j. Here each spin value 
is Sj = 1/2 and Sn = S . We restrict to the case of 
K 2 > — \K±\, so the ground state of H nec ki a ce is not 
ferromagnetic. We further specialize to N odd so the 
ground state is characterized by a spin quantum number 
S z = ±1/2. For K± < the total momentum of each of 
the S z = ±1/2 states is zero and the ground state is two- 
fold degenerate. For K ± > the ground state also has 
non-zero total momentum P = k\P\, k = ±1 implying 
four-fold degeneracy. In this paper we focus on the limit 
in which the Kondo coupling, J, of the cluster to the con- 
duction electrons is weak enough and the temperatures 
of interest low enough, that it is sufficient to consider 
only those spin flip processes that couple the degenerate 
ground states of the cluster. We expect that stronger 
couplings, implying mixing of higher energy states of the 
cluster, give rise to qualitatively similar physics. This 
case will be considered in a future paper [12]. 

If H^ upHng = 0, then the z— component of the cluster 
spin is a good quantum number. The cluster remains in 
a definite spin state, the conduction electrons just feel 
a potential due to the longitudinal coupling J z , and the 
Hamiltonian H = H neck i ace + H cond + H* oupHng is eas- 
ily diagonalizcd. For non-zero J± the problem is no 
longer trivially solvable. To study it, we use a pertur- 
bative expansion in Jj_ to establish a mapping onto a 
Coulomb gas, following the approach introduced by An- 
derson and Yuval [6]. Let denote the ground state 
of Hq when the cluster is in the S z = +1/2 state, |ft) and 
consider the ground-state-to-ground-state amplitude (al- 
ternatively the partition function in the imaginary time 
formulation) given in the interaction representation by 



Z(t) = * 



exp[i f H'(t')dt'] 
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where H'(t) = e- iHot H^ oupHng e iHot and T is the "time 
ordering operator" . Restricting to the ground state sub- 
space of the cluster, and expanding to 0(H' 2 ), 
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where Ry = Ri R J; t\ 2 = t\ — t 2 and 
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The expressions above are similar to the ones introduced 
in the context of the X-ray [13] [14] and the Kondo prob- 
lem [6], but with new features arising from the spatially 



extended nature of the magnetic object. The most im- 
portant of these new features are the space dependence 
of F, arising from processes in which an electron is cre- 
ated at one site and destroyed at another, and a nontriv- 
ial dependence of matrix elements such as (f|-| 5+ |J|) on 
necklace size and symmetry. We discuss these in turn. 

Following Rcf. [13] we express F as the product of 
open line part g = Fj (S) and a closed loop part (S) 
= e c which may be computed by an expansion in J z . 
The crucial object in this expansion is the Green's func- 
tions G(r, t) of H con d- We require the long time limit 
(t R/vf with R the radius of the necklace) in which 
G(r,t) —> pFG(r)/it : for a spherical Fermi surface 
G(r) = 2sa\{kFr)/vpr. In order to obtain a closed form 
for g we exploit the specific geometry of the necklace. 
We parametrize the positions of the spins by the az- 
imuthal angle <j>j = 2ir j/N, j G [0,7V — 1]. Using R^i = 
2i?sin(0 fe j/2), where (j) H = (j> k - <f> u G(R U ) G(<t> k i). 
We finally define the one-dimensional Fourier transform 
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where we choose m = — — ...0, ...^j^ to make the 
symmetry of the problem most apparent: G_ m = G m . 
Using the to— basis and defining V = —2J Z (ffl S z |ff), 
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where r c is the short-time cut-off separating two spin-flips 
and we have neglected terms of 0(J Z ) in the exponent. 

We now consider the matrix elements. We first note 
that if K± > then the necklace ground states |fr), |ff) 
have a nonvanishing momentum ±P$,i}. and the matrix el- 
ement depends on the momentum transfer: (-f|-| Sj~ |JJ.) = 
(ft | S+ |ft) cxp[iQj] where Q = (P^-P t ). When K± > 0, 
the absolute value of the matrix element decays faster 
with N when P t = Pp For K ± > when P t = -Pjj. 
and for K± < we find, defining a < 1, 



s+|ft>| =|(ftl^o + l^)h^ 
l^lft)| =l(ftl^|ft)|~^. 



(15) 
(16) 



For \K Z \ < \K±\, in terms of the Jordan- Wigner 
fermion representation [17], the states |ft) and |ft) have 
different boundary conditions, one periodic and the other 
antiperiodic. Considering the low energy form of the 
Hamiltonian given by the Luttingcr (massless) model 
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and adopting the "boundary-condition changing opera- 
tor" formalism [15] [16] we find 

a =i(l-i cos- 1 for\K x \<\K ± \ (17) 

which we verify by a detailed numerical analysis (see 
Fig. 2) of the exact lattice solution for the K z = case 
(xx— model [17]) for which a = 1/4. 
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Fig. 2 Plot of absolute value of M = Sj evaluated nu- 
merically for the K z = case as a function of N . A M is the 
matrix element when the Jordan- Wigncr fcrmion is added (when 
spin is flipped) in the unoccupied momentum fc M , each of which 
contributes to M. The inset plots its real, imaginary parts and the 
absolute value as a function of fe M for N = 283. The maximum 
contribution is near the Fermi point. We find that \M\ is equal to 
|Arj| (obtained when fc M is closest (fi = 0) to the Fermi point) and 
when plotted against N, the two plots overlap. 

For the Ising case (K z = oo) for K z > 0, the energy is 
minimized by having alternating up and down spins, but 
N being odd there is one frustrated bond with N possible 
locations for each of S z = —\ and |. A perturbative 
analysis around this 2N— fold degenerate point yields 



1 for JT Z » |J£i|. 



(18) 



We expect a = 1 for all K z > \K±\ (Ising regime). 

The results for a indicate that it is much easier than 
one might guess, to flip a large cluster between its two 
ground states. Quantum fluctuations, which for example 
in the Ising regime delocalize the frustrated bond, bring 
about dramatic enhancement in the tunnelling caused by 
one spin flip by a conduction electron. By distributing 
S z = ±1/2 over N sites, they also decrease the expecta- 
tion value of SJ, thereby reducing the potential felt by 
the conduction electrons at each spin site. 

We are now in a position to present our result for 
the ground-state-to-ground-state amplitude for the spin 
necklace embedded in a conduction bath, 
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where using 2irQ/N = Q = (Py. - P t ), 

£rn = Jz[G m + G m+ g] - J z ^ G 2 m , (21) 

m' 

u m = J±G m G m+ Q (22) 

with Jj_ = J±p F \(1l\S£ |JJ.)| - J ± p F /N a and T z = 
2p F J z (1t\ S§ ~ J zPF /N. Note that since J2 m G 2 m ~ 
N, the two terms in e m are of the same order in N and its 
leading behavior is given by the term linear in J z p F (as- 
sumed < 1). In the corresponding expression obtained 
in Ref. [11] it was assumed that only one channel was 
important and in the expression for e m , the term lin- 
ear in J z was overlooked. The above expression for Z, 
Eq.(19), is the second-order term in the expansion of a 
multicomponent Coulomb gas model in which each spin 
flip event (Coulomb gas charge) is labelled by the channel 
(to— value) of the electron- hole pair created. Higher order 
terms in the expansion in J± may be treated similarly, 
but compact expressions cannot be obtained because (as 
noted in a related context [18]) in the multi-channel case, 
the higher order terms do not combine into the Cauchy 
determinant form [6]. However we may define a scaling 
procedure as usual [19] by eliminating close-pair spin flips 
of the same channel. To the leading order the different 
channels are not coupled, and standard arguments show 
that near the weak coupling fixed point, the partition 
function retains the form Eq.(19) but with parameters 
E m (r c ) and u m (r c ) which evolve according to 



de m /d(\xiT c ) = 2u 2 m 
du m /d(\xiT c ) = u m e m /2 



(23) 
(24) 



with initial conditions given by Eq.(21) and (22). The 
flow implied by Eq.(23) and (24) is shown in Fig. 3. 




Fig. 3 Plot showing the flow for each channel, near the weak 
coupling fixed point. The heavy solid line is the separatrix e m = 
—u m /2. The shaded oval region shows the generic initial conditions 
for the xx symmetry. The dotted circle shows the generic initial 
conditions for the Ising limit. The striped region near u m = is 
the range of initial conditions for 'imperfect' necklaces. 

For each channel the flow is to the no-flip fixed line 
(u m = 0) if s m < —u m /2 and towards the strongly 
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coupled (rapidly flipping) regime (with e m — u m /2) if 
£m > —Mm/2- Scaling amplifies the differences between 
channels. For example, in the K± < case of the 
xx— model, the initial condition is tt TO = J ± ppGm/N 1 / 4 
and £^ = 2J z ppG m /N. For large N, is thus para- 
metrically smaller than u^; all initial conditions are near 
the vertical axis and the dominant channels are those 
with the largest G m . Representative results are shown 
in Fig. 4. Note that for the ideal ring studied here, 
G- m = G m . We see that all channels with m < Nkp/^ir 
couple with roughly equal strength, but generically there 
is one largest pair, so after a transient many channel 
regime, the asymptotic behavior is of the two-channel 
Kondo problem. Exceptions can occur, and in particu- 
lar for fcp < 2it/N, the m — channel dominates. The 
same generic behavior occurs in other cases we have stud- 
ied, but we note that in the Ising regime the J z and J± 
couplings are both ~ 1/N, so the initial conditions are 
shifted as indicated in Fig. 3. 
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Fig. 4 Plot showing G m in the K± < case of the xx— model 
for two values of kp . For kp = w, the pair m = — 10 and 10 and 
for kp = 7r/2, the pair m = —4 and 4 is the highest, constituting 
the dominant pair which governs the low energy behavior. 

To summarize: we have shown that the generic long- 
time dynamics of a spin necklace coupled to a conduction 
electron bath is of the two-channel Kondo class. We em- 
phasize that generic deviations from perfect translation 
invariance break the two-channel symmetry. A key fea- 
ture of the model is the large quantum fluctuations char- 
acteristic of the isolated necklace: these ensure that the 
dominant spin flip amplitude between the ground states 
of the spin necklace decays only as a power of the sys- 
tem size (N~ a with < a < 1) and that the static 
field exerted by the necklace on the conduction bath is 
small, of order 1/N, so "orthogonality" effects do not 
dominate. To explore the changes occurring as quantum 
fluctuations are reduced, we consider a necklace in the 
extreme Ising limit (K z 3> and with one weak link 

(K® 1 — > K z (l — 5)) which acts to localize the frustrated 
bond. A detailed discussion and extension to other cases 
will be published elsewhere [12]. Here we note that al- 
though the residual reflection symmetry allows two chan- 



nel behavior, the spin flip amplitude changes from ~ 1 /N 
to ~ exp[—N]n(K±/6K z )] (moving the initial condition 
much closer to the u m — axis of Fig. 3). More im- 
portantly, because the frustrated bond is localized, the 
quadratic, "orthogonality" term in e m (Eq.(21)) becomes 
~ J Z N, moving the initial condition substantially to the 
left (independent of the sign of J z ); and, for large enough 
N, into the stable (no asymptotic flip) region. 

Extensions of the present work include a treatment of 
systems with open boundary conditions, and of other 
nanostructures. We expect generically that internal 
quantum fluctuations of a nanoscale object (neglected in 
many treatments, Ref. [8] and [11] for example) will have 
a dramatic effect on the long time dynamics. 
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